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The remarkable phenomenon of weak magnetization hysteresis loops, observed recently deep in 
the vortex-liquid state of a nearly two-dimensional (2D) superconductor at low temperatures, is 
shown to reflect the existence of an unusual vortex-liquid state, consisting of collectively pinned 
crystallites of easily sliding vortex chains. 



Many potentially important superconductors, such as 
some of the high-T c cuprates, as well as the ET-based or- 
ganic conductors, are extremely type-II superconductors 
with small in-plane coherence lengths and nearly 2D elec- 
tronic structure. Thereby, ET (or BEDT-TTF) stands 
for bisethylcncdithio-tctrathiafulvalene. Drastic devia- 
tions from the predictions of the mean-field theory for 
these materials, due to strong fluctuations in the super- 
conducting (SC) order parameter, are therefore expected, 
in particular under strong perpendicular magnetic field 
Q. The great interest, from a fundamental point of view, 
in the latter class of materials stems from their relatively 
low upper critical fields, which enable to investigate the 
virtually unexplored phase diagram and vortex dynamics 
of strongly type-II superconductors at low temperatures 
and high magnetic fields. 

Of special interest is a recent striking observation in 
/3"-(ET) 2 SF 5 CH 2 CF 2 S03 of small, but very clear magne- 
tization hysteresis loops appearing well above the 'major' 
irreversibility field at low temperature, where significant 
de Haas-van Alphen (dHvA) oscillations are observable 
as well @ (Fig. QJ. It should be stressed that the oc- 
curence of such high field hysteresis tail is not peculiar 
to this particular material. It can also be observed at low 
temperatures in k-(ET) 2 Cu(SCN) 2 0. 

The unusual feature of this irreversibility effect is as- 
sociated with its appearance deep in the vortex-liquid 
phase, where one usually expects unrestricted motion of 
flux lines through the entire SC sample, ensuring the 
establishment of thermodynamic equilibrium. In the 
present letter we argue that the observed hysteresis ef- 
fect is a general feature which reflects the unusual nature 
(namely the nematic liquid-crystalline structure) of the 
low-temperature vortex-liquid state in 2D superconduc- 
tors well above the vortex-lattice melting point, which 
was recently predicted theoretically . 

According to this model the vortex state above the ma- 
jor irreversibility field, Hi rr , (which is found to be close 
to the vortex-lattice melting field, see below) is not an 
isotropic liquid but some mixed phase containing SC do- 
mains of easily sliding parallel vortex chains, which are 
stabilized by a small number of strong pinning centers. 



This model resembles the sluggish vortex fluid picture, 
proposed by Worthington et al to describe the inter- 
mediate vortex liquid phase in defect-enhanced YBCO 
singel crystals, and applied more recently by Sasaki et 
al [f| to account for a similar behavior in the organic 
superconductor k-(ET) 2 Cu(SCN) 2 . As shown below, a 
simple Bean-like model for the magnetic-induction pro- 
files, associated with the injected vortex chains into the 
SC sample, yields very good quantitative agreement with 
the measured field dependence of the magnetization hys- 
teresis. 

Our model consists of independent 2D SC layers in 
the x-y plane under the influence of an external uniform 
magnetic field H = HH, H > 0, perpendicular to the 
layers (Fig. |2J. The SC sample is in contact with a nor- 
mal metal at two parallel planes, x = 0, and x = L x . 
H is varied, at low temperatures, from above H C 2 to the 
irreversibility field Hi rr < i? c2 , where magnetic quan- 
tum oscillations can be detected. The dHvA effect is 
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FIG. 1: Field dependence of the torque signal with several 
small hysteresis loops in the tail of the hysteresis loop shown 
in the inset. Arrows indicate the field-sweep directions. Inset: 
Lower field torque data showing how the major irreversibility 
field is determined. 
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measured by means of the torque method, in which the 
signal is detected during steady (upward and downward) 
sweeps of the external magnetic field. Note that, for the 
sake of simplicity, the small in-plane component of the 
external magnetic field, required for this high-resolution 
measurement, is neglected in our analysis. 

The influence of a steady increase (or decrease) of the 
magnetic field on the distribution of magnetic flux lines 
within the SC region is determined by the pinning forces 
acting on quantized SC vortex lines near the normal- 
superconducting (N-S) boundary planes. The pinning- 
force resistance against flux injection leads to the es- 
tablishment of a flux-density gradient perpendicular to 
the magnetic field direction along the normal of the N-S 
boundary plane (along x). Assuming the flux injection 
to be uniform along y, the current density j = jy and 
the field gradient, are connected by Ampere's law: 
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where the Lorentz- force density (per unit volume), ex- 
erted on the vortex current by the magnetic field, is 
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For magnetic fields below Hi rr , 
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the pinning force, 
pm — , is sufficiently strong to balance the ac- 
tion of the driving Lorentz force, Fl = \jB = — 
with U L (x) = ±[B(x)} 2 , so that ^f- = 0, with 
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Utot{x) = Upi n (x) + Ul{x). Consequently, the static 
magnetic- field gradient, established across the sam- 
ple, is followed by a dissipationless current parallel to the 
surface. This static equilibrium situation persists under 
the increasing local magnetic potential Ul(x) near the 
surface, up to the 'elastic limit' of the pinning poten- 
tial, Upi n (x), above which bundles of vortices, accumu- 
lated near the N-S boundary, can jump over (or tunnel 
through) the £/t t-potential barriers to neighboring po- 
tential wells. Such collective vortex jumps into the SC 
region quickly reduce the local Lorentz force by reliev- 
ing the large local field gradient, so that the movement 
stops after a few steps. The observed sharp magnetiza- 
tion jumps in this field range 0i0 seem to be associated 
with this effect. 

The vortex-vortex correlations in this field range are 
relatively strong and long ranged, so that the effective 
pinning force is long ranged as well. The effect is par- 
ticularly significant at low pinning-center density, where 
these correlations can dramatically enhance collective- 
pinning processes, since many vortices are involved in the 
interaction with a single pinning center. Furthermore, 
due to the strong anisotropy of the vortex-vortex cou- 
pling Q, this collective-pinning mechanism is also very 
anisotropic. Thus, the magnetic- flux injection should en- 
counter different resistances at different locations of the 
N-S boundary, depending on the local orientation of the 
vortex crystallite with respect to the surface. The most 
favorable orientation is that of crystallites having their 
principal lattice vector aligned perpendicular to the N-S 
boundary, i.e., along x, where a weak driving force can 
inject many vortex chains into the SC region. 

As the external magnetic field approaches Hi rr from 
below and the strength of the pinning potential is de- 
creasing (because U P i n is proportional to the local su- 
perfluid density, which decreases with increasing B), the 
increasing Lorentz force can exceed the maximal value of 
the pinning force, i.e., the highest maximum of Ut t(x) 
turns into an inflection point. Above this critical field, 
Utot (x) has no barriers and the motion of vortices in the 
entire SC region becomes a continuous flow. This vortex 
movement with velocity is opposed by a dynamic fric- 
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of the driving Lorentz force at the critical velocity v% 
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In the corresponding steady state the motion of flux lines 
generates an electric field 



FIG. 2: (a) Sketch of a vortex crystallite, oriented with its eas- 
ily sliding Bragg chains perpendicular to the N-S boundary, 
under injection of flux lines. Note the scaling of the magnetic 
length over a macroscopic distance Lq. (b) Magnetization 
profiles during upward and downward field sweeps. 
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with E = \Bv c ^ parallel to the current density j. Con- 
sequently, the system develops a finite electrical resis- 
tivity to the critical-current flow, j c = crjv^/B, namely 
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p[B) = E/j = B 2 /c 2 r]. Now, the Bardeen-Stephen (BS) 
relation yjll0|, V * BH C 2/c 2 p n , which is equivalent to 
the simple linear magnetoresistivity || p = (B/H C 2)p n , 
together with Eqs. Q and Q, lead to: 
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The Faraday induction law together with Eq. yield 



the flux-line conservation law, 



OB 

at 



0. which 



may be combined with Eq. (J5J to obtain the nonlinear 
'diffusion' equation for the magnetic-flux density: 
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with the diffusion coefficient = ^f-. In the high- 
field limit of interest here, when the range AH of a sin- 
gle field-sweep cycle is much smaller than the initial field 
of the hysteresis loop, Hq, the nonuniform part of the 
magnetic induction associated with the flux flow is very 
small, i.e., \b(x,t)\ = \B(x,t) - H(t)\ < H{t) ~ H . Un- 
der this condition the second term on the RHS in Eq. © 
can be neglected, and the equation may be linearized to 

W ~ (:lcir) D <l>7%r- The solution of tnis equation, sat- 
isfying steady upward and downward field-sweep bound- 
ary conditions, ^{x = 0,t) = ^-{x = L x ,t) = 
for < t < t, and for r < t < 2t, takes the 
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forms: B + (x,t) = Hq 
B_(x,t) = H + AH (2 
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spectively, with r half the period of the field-sweep cycle. 

Note the spatial rigidity of the induction profiles 
B±(x,t), which is assured, in the high-field limit 
\b(x, t)\ <C H, by the nearly instantaneous propaga- 
tion of any local fluctuation of the magnetic-flux den- 
sity. Indeed, the propagation velocity of such a fluctu- 
ation, Vf = 7jf/§f , is found by Eq. 10 to be: Vf = 



-— ( B^r) /H , so that the flux- line velocity ty, (see 
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Eq. ©) satisfies: \ V<f> \ ~ (ff) 2 /^ 

The resulting Bean-like induction profile can now be 
exploited to evaluate, for each small hysteresis loop, 
the jump, AMfj,, of the magnetization occurring at 
a point where the field sweep is reversed (Fig. 0). 
The corresponding jump, AM^j. = -^AB, is deter- 
mined by the spatial average of the induction differ- 
ence, [B_ (x,t) — B + (x,t)], where the (field-sweep rate- 
dependent) parameter £ (( < I) is an effective volume 
fraction of the SC domains in the sample. This simple 
Bean-like model may be tested by estimating the parame- 
ter p„, which should be identified with the in-plane resis- 
tivity at H C 2 of the 2D conductor under study. This can 
be done by using the experimental values: AirAM^i sa 4x 
I0" 7 T (sec below), L x w 0.1 cm, and (A£) « 0.005 T/s 
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FIG. 3: Experimental major irreversibility fields (open cir- 
cles), and theoretical melting fields (solid lines), as functions 
of temperature for two organic superconductors. Data for 
k-(ET) 2 Cu(NCS) 2 from Ref. 0. 



in the relation {AH/t)( = BD^Hq/ H c2 )AB / L 2 X , 



which yields p n 
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3c 2 AM u V t , 

Assuming that C, is of the order unity, this is a reasonable 
result, (i.e., of the order of a typical metallic resistivity), 
which cannot be easily verified experimentally, since for 
the very anisotropic metal under study here, measure- 
ments of the in-plane resistivity are extremely sensitive 
to out-of-plane scattering by a minute amount of defects. 

It is interesting to estimate the characteristic energy 
scale corresponding to the friction of moving vortices 
in the liquid state above Hi rr against the effective pin- 
ning force. This can be done by estimating the total 
spatial variation of the magnetic induction within the 
SC sample, £;c|§§| ~ CAB, with the bar indicating 
spatial averaging, and comparing it to the experimen- 
tal value of AMjj. Thus, the work (per unit volume) 
done by the driving Lorentz force against the friction 
force, E frlc = \ 3c BL x = &B\%f\L x « ^BAB, may 
be rewritten in units of the maximal SC condensation- 
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energy density, E cond 

s fric = E fric /E cond ?s 4k 2 (AB/H c2 . Using the experi- 
mental value of the magnetization jump at 2 T, we find 
(AB = 37rAA% w 4x 10" 7 T, so that e fnc w 1.7 x 10" 3 , 
where k « 46 is used 01 • This estimate is smaller then 
the characteristic energy associated with the melting of 
the vortex lattice in 2D superconductors, which is about 
2% of the SC condensation energy. It may indicate that 
most of the vortex crystallites are aligned with their prin- 
cipal lattice vector along the direction of the flux-line 
motion so that the flux flow is dominated by eas- 
ily sliding vortex chains, which are subject only to the 
residual shear resistance characterizing the vortex-liquid 
state above the melting transition This picture is 
very similar to the moving smcctic state found recently 
in many numerical simulations of 2D vortex lattices in 



4 




B(T) 



FIG. 4: Field dependences of the experimental magnetization 
jumps (triangles) and the fit (solid line) according to Eq. J7J. 



the presence of random pinning centers |l3j |. 

The value of Sf r ic estimated above may also reflect 
the correlation between the onset of irreversibility and 
the vortex-liquid freezing transition. This may be ver- 
ified by calculating the temperature-dependent melting 
field, H m (T), and comparing it to the irreversible field, 
Hi rr (T), extracted from the experimental data. Similar 
to Ref . |3j , Hi rr (T) is determined here by the onset of the 
major hysteresis in the magnetization curve ( see Fig. ^ 
). Within the Ginzburg-Landau (GL) approach, H m (T) 
has been derived by several authors [4J, |l4j . It can be ob- 
tained from £ 2 (iJ m ,T) = g 2 n , where £ 2 (H,T) = e /k B T, 
with £q = a 2 /2f3 being the SC condensation energy per 
unit vortex, and a and (3 are the GL parameters 
Here we estimate |A | 2 = f = (1.76k B T c ) 2 In (^), and 

P = wtj^JZp ■ Usin S our estimate, g 2 m w 1/4A 2 , with 
A ~ 0.066 the effective minimal shear energy relative 
to £0 U and neglecting the weak temperature depen- 
dence of H C 2 at low T, we obtain the simple equation 
for h m (t) = H m (T) I H c2 (0): ln(/i rn ) = -k Q Vth m , where 
t = T/T c (0) and fco is a single dimensionless parame- 
ter depending on the properties of the superconductor 
through k 2 ~ Q.Vbg 2 n E F (Tiuj C 2) 2 / \k B T c f . The function 
h m (t) for /3"-(ET) 2 SF 5 CH 2 CF 2 S0 3 , with T c (0) = 4.4 K 
and H C 2 = 2 T is shown in Fig. OJa). The experimen- 
tal irreversibility line agrees rather well with the calcu- 
lated melting curve. A similar procedure is applied to 
k-(ET) 2 Cu(SCN) 2 . The parameter H c2 w 4.7 T was de- 
termined by fitting the additional damping amplitude of 
the dHvA oscillation in the mixed SC state, as calculated 
from our SC fluctuations theory 0, to the correspond- 
ing experimental data 0. Again, the calculated melting 
curve is found close to the irreversible line obtained in 
Ref. [Fig.EJb)]. It is interesting to note that at T = 
the resulting melting field coincides with the mean-field 



H C 2- The sharp increase toward -ff c2 , characterizing H m 
at T — > 0, is similar to that of Hi rr observed experimen- 
tally in both materials. 

The relation H irr ~ H m motivates us to exploit our 
fluctuating-vortex-chain model above H m in calculating 
the magnetic-field dependence of the jump AM^ above 
H irr ) and compare it to the experimental result (Fig.nj. 
Combining Eq. with Eq. © , we may rewrite the crit- 
ical gradient as \dB/dx\ = AttJe^/B = 4irE fric /L x B. 
Since the generators of the friction force, i.e., the pin- 
ning force and the vortex-vortex interactions, all origi- 
nate in the existence of local SC order at the vortex posi- 
tion, the friction energy Ef r i C should be proportional to 
the mean-square SC order parameter, (|A (-ff)| 2 ), and 
so the magnetization jump can be written as AM^ = 

±AB » ^j|f[ cx ±(\A (B)\ 2 ). Using the limit of 



independently fluctuating vortex chains in the GL the- 
ory [l| to describe the vortex- liquid state above Hi rr , 
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is a universal function of the dimensionless parameter 
£(-B, T). Employing the material parameters, T c = 4.4 K, 
— = 2.0, Ep/ks = 133 K at the temperature of the ex- 
periment, T = 40 mK, so that £ ks 31(1 - B/H c2 )/B, 
and treating if c2 and the proportionality factor in Eq. 
(J7J as adjustable parameters, the best fit is obtained 
for H C 2 = 2 T (Fig. 0J. The resulting curve reflects 
the crossover between the vortex-crystal state below H c2 , 
which is well described by mean-field theory, and the nor- 
mal state far above -ff c2 , with the long tail of the field- 
dependent magnetization hysteresis corresponding to the 
enhanced influence of 2D SC fluctuations. 

In conclusion, we have shown that the striking phe- 
nomenon of the weak magnetization hysteresis loops, ob- 
served deep in the vortex-liquid state of a nearly 2D su- 
perconductor can be reasonably explained as arising from 
shear viscous flow of easily sliding vortex chains, which 
are clustered around a small number of strong pinning 
centers. 
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